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Recent developments in finite-temperature studies of lattice QCD are reviewed. Topics include (i) tests of 
improved actions for the pure gauge system, (ii) scaling study of the two-flavor chiral transition and restoration 
of Ua{1) symmetry with the Kogut-Susskind quark action, (iii) present understanding of the finite-temperature 
phase structure for the Wilson quark action. New results for finite-density QCD are briefly discussed. 



1. Introduction 

Finite-temperature studies of lattice QCD have 
been pursued over a number of years. Quite 
clearly the pure gauge system is the best under- 
stood of the entire subject. The system has a 
well-established first-order deconfinement transi- 
tion[0|, and extensive and detailed results are al- 
ready available for a number of thermodynamic 
quantities!^. Nonetheless many new studies have 
been made for this system recently. The purpose 
is to examine to what extent cutoff effects in ther- 
modynamic quantities are reduced for improved 
actions as compared to the plaquette action which 
had been used almost exclusively in the past. 

Full QCD thermodynamics with the Kogut- 
Susskind quark action has also been investigated 
extensively in the past. A basic question for this 
system is the order of chiral phase transition for 
light quarks. For the system with two flavors, 
finite-size analyses carried out around 1989-1990 
indicated an absence of phase transition down to 
the quark mass rUq/T « 0.05^, and a more re- 
cent study[Q| attempted to find direct evidence for 
the second-order nature of the transition, as sug- 
gested by the sigma model analysis in the contin- 
uum[^, through scahng analyses. Scaling studies 
have been continued this year to establish the uni- 
versality nature of the transition on a firm basis. 
Another issue discussed at the Symposium is the 
question of restoration of Ua{^) symmetry at the 
chiral transition. Results have been presented for 
equation of state both without and with use of 
improved actions. 

Studies of thermodynamics with the Wilson 



quark action is much less developed compared to 
that for the Kogut-Susskind quark action. Past 
simulations found a number of unexpected fea- 
tures, which made even an understanding of the 
phase structure a non-trivial problem|^. Re- 
cently, however, considerable light has been shed 
on this problem through an analysis based on the 
view that the critical line of vanishing pion mass 
marks the point of a second-order phase transi- 
tion which spontaneously breaks parity and flavor 
symmetry [0. Some new work with improved ac- 
tions, which was initiated a few years[^, has also 
been made this year. 

In this article we review recent studies of finite- 
temperature lattice QCD. In Sec. 2 we summarize 
results for the pure gauge deconfinement tran- 
sition obtained with a variety of improved ac- 
tions. Results for the two-flavor chiral transition 
for the Kogut-Susskind quark action are discussed 
in Sec. 3 with the main part devoted to scaling 
analyses of the order of the transition. In Sec 4 
we describe recent progress on the phase struc- 
ture analysis for full QCD with the Wilson quark 
action. This year's results for finite density are 
briefly discussed in Sec 5. Our summary and con- 
clusions are presented in Sec. 6. 

2. Recent vifork on pure gauge system 

In Table |l| we list recent studies of the pure 
gauge system using improved actions. Among ac- 
tions constructed through renormalization group, 
RG(1,2)[Q includes 1x2 loop in addition to the 
plaquette. The action FP is an 8-parameter ap- 
proximation to the fixed point actional contain- 
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Table 1 

Recent work on pure gauge system with improved 
actions. Argument for [3c means spatial volume, 
n{L) the torelon mass for spatial size L and Nt 
the temporal lattice size. 



action 



ref. 



measurements 



Nt 



RG-improved 

FP (type I) |L 
FP (type Ilia) |10| 
FP (type Ilia) illj 
RG(1,2) |12| 
|13| 



Symanzik- improved 



^c(oo),/i(i) 

/3c(oo) 

p 

l3c{oo),a 
/3c(oo),/i(i) 



2,3,4,6 
2,3,4,6 

2,3 

4,6 

2,3 



S(l,2) 



tree 



|14| 
|15| 

S(l,2)tadpole Hi 
S(2,2)^;j.ge 

SLWtadpole 



Pcim) 3,4,5,6 
Pc{co),a,e,p,ai 4 
/3c(oo),cr, e,p, cr/ 4 
[i6f,||/3,(oo),a,6,p 4 
[17| (3c{2Nt),^i{L) 2,3,4 



13| Pci^),KL) 



2,3 



ing plaquette and unit parallelogram and their 
power up to 4; type I |^ and type llla jl^ differ 
in the blocking procedure, the latter including 5 
and 7 link staples to improve rotational symme- 
try. Actions improved according to Symanzik's 
program with the addition of an m x n loop are 
denoted as S(m,n), distinguishing tree- level coef- 
ficients from tadpole-improved ones ||l9[| by corre- 
sponding sufhces. The S(l,2) action improved to 
one loop order ll^] including tadpole factors is 
denoted by SLWtadpole- 

2.1. Critical temperature 

2.1.1. r,/v^ 

A basic quantity for the pure gauge system is 
the ratio Tc/ \fa ^ where cr is the string tension ex- 
tracted from the static quark-antiquark potential. 
In Fig. |l| we plot recent results for this ratio for 
improved actions and the plaquette action]^ as 
a function of (aTc)^ = l/A^t^ with Nt the tempo- 
ral lattice size. Results are not available for fixed 
point actions. For SLW^^dpole' ^'^'^ which values 
of fj are not given by the authors [pT| , we use an 
interpolation of the results reported in ref. [^^ . 

Scaling behavior can be examined for two types 
of improved actions, SLWtadpole ^-i^d RG(1,2). 
Comparing their results with those of the pla- 
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Figure 1. for pure gauge system normalized 
by square root of string tension -^/ct as a function 
of (aTc)^- Dotted fines are 0(c?) fits. 

quette action, it is apparent that the improved 
actions exhibit a better scaling behavior. For 
SLWtadpole ratio increases by about 5% over 
aTc « 0.25— 0.5, while a variation of 3— 7% occurs 
for the plaquette action over a factor two smaller 
range aT^ ~ 0.08 - 0.25. The resuhs for RG(1,2) 
are constant over aTc ~ 0.17 — 0.25 within the 
quoted error of 1 — 2%. The errors are still sizable, 
however, to draw a conclusion on the magnitude 
of slope from the limited range of lattice spacing 
explored so far. 

We observe, however, that the continuum ex- 
trapolation of Tel \fo for the improved actions 
do not agree with that for the plaquette action. 
For SLW^adpolei assuming an O(a^) or O(a^) 
dependence (we ignore the factor of the ac- 
tual form 0((7'*a^)) both of which are consistent 
with present data, we find Td \fo — 0.641(7) and 
0.653(10) in the continuum limit. For RG(1,2) 
we expect an Oic?) scaling violation, with which 
we obtain Td \fd = 0.663(13). These values are 
2 — 5% (one to two standard deviations) larger 
compared to the estimate for the plaquette ac- 
tion Tc|^Jo = 0.629(3)@ (open circle at aT^ = 0) 
obtained with a quadratic extrapolation. 

We stress that efforts to resolve the discrepancy 
should be made for we would then have a deter- 
mination of the basic ratio T^j \fo accurate at the 
level of 1 - 2%. 

Let us add a remark on the results at 
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oTc — 0.25 where values for six types of ac- 
tions are available. Among those belonging to 
the category of Symanzik improvement, we ob- 
serve that Tcl\fo systematically increases in the 
order, plaquette^S(l,2)tree ^S(l,2)tadpole ^ 
SLW^j^jjpojg, reaching T^j \/a sa 0.64. If we take 
Tel ~ 0.65 as the continuum value(see above), 
this trend is consistent with the theoretical expec- 
tation that cutoff effects are reduced with an in- 
creasing degree of improvement. In this regard 
S(2,2)tree seems less improved than S(l,2)ti-ee- 
The fact that the value for RG(1,2) lies above 
that for SLW^j^fjpQig may be ascribed to a larger 
coefficient of the 1x2 loop term for this action 
compared to that for the latter. 

2.1.2. Te/v/^ 

Another quantity often used for testing im- 
provement with simulations on small lattices 
is the torelon mass extracted from the 

Polyakov loop correlator on a lattice of spatial size 
L. Defining ct(L) = ^{L)/L we compile in Fig. ^ 
results for the ratio Td \J <j(V) for a fixed physi- 
cal spatial size L = 2/Tc as a function of (aTc)^. 
The results for the fixed point action FP(type I) 
and the plaquette action were already available 
last ye arpd except for the value for the former for 
TVt = 6|23|]. This year, results have been reported 
for the actions SLW^adpole ^"^^ RG(1,2). 

Taking all the data together, one finds a bet- 
ter scaling behavior exhibited by the improved 
actions compared to the plaquette action. On 
a closer look, however, some discrepancy is ob- 
served among data for improved actions. While 
results for FP(type I)Q show a constant behavior 
in aTc within the quoted error of 1 — 2%, those 
for SLW^adpole obtained by the Cornell group [0 
exhibits an increase of 4% over aT^ = 0.5 — 0.25, 
and increasingly deviate from those of FP (type 
I) toward smaller lattice spacings. The results for 
RG(1,2) lie in between those of the two actions. 
There is also a discrepancy among results from 
ref. and II, both for SLWtadpole- 

In order to extract the torelon mass, the Cor- 
nell group [r^ employs a multi-state fit including 
excited states to a set of Polyakov loop correlators 
unsmeared or smeared at either source or sink, 
while a single state fit to a correlator smeared at 
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Figure 2. normalized by \J o{V) as a function 
of (aTc)^ where a[L) = n{L)/L with yLt(L) the 
torelon mass for a physical spatial size L fixed at 
L = 2/Tc. Dotted lines are 0{a?) fits. 

source and sink is used in the other studies |^,|l3|. 
The discrepancy possibly originates from contam- 
ination of excited states in the latter fit||l7[|. 

Because of the different trend in the lat- 
tice spacing dependence, continuum extrapola- 
tion leads to a significant scatter of the ratio 
among the actions. Dotted lines in Fig. ^ il- 
lustrate the difference, where we assume O(a^) 
dependence since O(a^) terms are expected to 
be present for all the actions employed. Com- 
pared to the value for the plaquette action, the 
values for SLW^^dpole RG(1,2) deviate at a 
±(1 - 2)% level, while that for FP(type I) is 6% 
lower. 

The discrepancy should be resolved, particu- 
larly to ascertain if cutoff effects are reduced 
to within 1 — 2% at a large lattice spacing of 
aTc — 0.25 — 0.5 for the fixed point action. 

Another problem concerns the consistency be- 
tween the values of T^j \Ja[V) and Td^fo. The 
relation a ~ ^{L)^^/ (3L^) which holds in string 
models [p4| suggests that T^j \foiJj) increases by 
about 0.04 compared to Tc/y/a for L = 2/Tc. 
Making a comparison for each action, we find that 
the increase for SLW^^dpole consistent with this 
estimate. For the plaquette action the value of 
Tcj \J a(L) in the continuum limit is 3% higher 
than the estimate, while for RG(1,2) the value is 
3% lower. 
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Figure 3. (a) latent heat Ae/T^and (b) interface 
tension <Ji/T^ for plaquette [p5| and Symanzik- 
improved S(l,2) action [p^. 

A systematic comparative analysis of various 
actions employing the same simulation and anal- 
ysis procedures is need to resolve this discrepancy. 

2.2. Latent heat and interface tension 

The latent heat Ae and interface tension cr/ at 
the deconfinement transition exhibit large scaling 
violations from an A^t = 4 to an Nt — Q lattice 
for the plaquette action as shown by open circles 
in Fig. ^. The Bielefeld group measured these 
quantities for the S(l,2) action||l5| both without 
and with tadpole improvement (triangles). 

For both quantities the rapid decrease of values 
for the plaquette action indicates that the contin- 
uum value would be lower than that on an A^t = 6 
lattice. The values for the improved actions are 
indeed smaller already for A^i = 4. 

2.3. Energy density and pressure 

Bulk thermodynamic quantities such as energy 
density and pressure receive substantial contribu- 
tions from high momentum modes. Since cutoff 
effects in these fluctuations are directly reduced 
by the improvement procedure, we may expect 
significant improvement. 

Numerical tests have been made[^,^ for the 
actions S(l,2)tj.ee, tadpole and S(2,2)tree through 
comparison of pressure calculated on an A't = 4 
lattice with an estimate of the continuum value 
obtained with the plaquette action^. The re- 
sults for the improved actions are close to the 
continuum estimate already for Nt — 4, especially 



for S(l,2)^g^(jpQjg for which a good agreement is 
seen even close to Tc- The fixed point action 
FP(type Ilia) also exhibits a similar agreement 
for A^t = 3[|ll- We refer to Figure 1 of ref. 
for these points. 

2.4. Summary 

Improved actions lead to a sizable reduction of 
cutoff effects in thermodynamic quantities. De- 
tailed scaling studies, however, are limited to that 
of critical temperature, for which consistency of 
results for various actions are not yet attained be- 
yond the level of 5%. Further studies are needed 
to see if accurate determination of thermody- 
namic quantities at a few percent level is possible 
through simulations with a moderately large tem- 
poral size of 4 < A^t < 8 with improved actions. 

3. Chiral phase transition with the Kogut- 
Susskind quark action 

3.1. Order of transition for Nf = 2 

In Table || we list major studies on the order 
of two-flavor chiral phase transition carried out 
on an A^t = 4 lattice with the spatial size Li^ . In 
the finite-size scaling study pursued around 1989- 
1990 [|2^ , ^ , measurements were made of the sus- 
ceptibility of the Polyakov line Vl given by 



Xn = [((Ref^)') - (Ref])2] 



(1) 



and a pseudo-susceptibility of chiral order param- 
eter defined by 



Xc 



1 

9L3 



(2) 



where ^ is a gaussian noise and D the Kogut- 
Susskind quark operator. It was found that the 
peak height of the two susceptibilities increases 

Table 2 

Studies of order of Nf — 2 chiral transition on an 
A^t = 4 lattice. 



ref. 



size 



KEK(1990) 

Columbia(1990) 

Bielefeld(1994) 

Bielefeld(1996) 

JLQCD(1996) 



[26| 6^,83,123 

[281 
[291 



16^ 

83 

123, 16^ 
83,123, 16^ 



0.0125,0.025 
0.01,0.025 
0.02- 0.075 
0.02- 0.075 
0.01 - 0.075 
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up to 12"^, but stays constant within errors be- 
tween 12'^ and 16'^ both at — 0.025 and 
0.0125 - 0.01 (see Fig. 11 and 12 in the second 
paper of ref. The saturation imphes the ab- 

sence of a phase transition down to niq « 0.01. 
Since this quark mass is quite small, correspond- 
ing to rriTr/mp w 0.2 at the point of the transition 
(3c ~ 5.27, it was thought that the result is con- 
sistent with the transition being of second order 
at = as suggested by the sigma model anal- 
ysis!]. 

3.1.1. Scaling analysis of susceptibilities 

One can attempt to examine if the transition is 
of second order employing the method of scaling 
analysis. Let us define the susceptibilities Xm and 

Xt,i {i = f, cr, t) by 

Xm = V[{iqqr)-{qqr] (3) 
XtJ = V[{qq-qD„q)-{qq){qDoq)] (4) 
XM = V[{qq-P,)-{qq){m, i = <t,t (5) 

with V — L'^Nt, Dq the temporal component of 
the Dirac operator, and Pg-^r the spatial and tem- 
poral plaquette. For a given quark mass m^, let 
(7~^(mq) be the peak position of Xm as a func- 
tion of the couphng constant g~^, and let Xm"^ 
and xTf^i^ — Z''^;''') be the peak height. For a 
second-order transition, these quantities are ex- 
pected to scale toward vriq — > as 



9c ^i^q) 



max 



max 
KtA 



f, <^,T 



(6) 
(7) 
(8) 



Let us note that xt,i (* = /, f, t) are three parts of 
the susceptibility Xt = V [(99 • e) — {q<i){^)] with 
e the energy density Q. The leading exponent is 
therefore given by Zt = Max{ztj, zt^a, -^t.r)- 

Natural values to expect for the exponents 
Zg, Zm and Zt at a finite lattice spacing are those 
of 0(2) K, U{1) corresponding to the exact sym- 
metry group of the Kogut-Susskind action. How- 
ever, sufficiently close to the continuum limit 
where flavor breaking effects are expected to dis- 
appear, they may take the values for 0(4) « 
SU{2) (g) SU{2) which is the group of chiral sym- 
metry for Nf = 2 in the continuum. One should 
also remember that mean-field exponents control 



the scaling behavior not too close to the transi- 
tion. A possibility of mean-field exponents arbi- 
trarily close to the critical point has also been 
discussed Q. 

The initial scaling study was carried out by 
Karsch and LaermannjJ] employing an 8^ x 4 lat- 
tice and niq = 0.02,0.0375,0.075. Compared to 
the 0(4) values their results for exponents show 
a good agreement of Zm, a 50% larger value for Zg 
and a value twice larger for zt- Comparison with 
0(2) and mean- field exponents is similar since 
they are not too different from the 0(4) values. 

This work had limitations in several respects: 
(i) the scaling formulae are valid for a spatial size 
large enough compared to the correlation length. 



At 



0.02 the pion correlation length equals 



sa 3. Whether the spatial size of L = 8 em- 
ployed is sufficiently large has to be examined, 
(ii) The size of the scaling region in terms of 
quark mass is a priori not known. Hence the 
behavior for smaller quark masses should be ex- 
plored to check if the results are not affected by 
sub-leading and analytic terms in an expansion 
of susceptibilities in ruq . (iii) In the original work 
the noisy estimator with a single noise vector was 
employed to estimate disconnected double quark 
loop contributions. This introduces contamina- 
tion from connected diagrams and local contact 
terms, which has to be removed. Other fac- 
tors such as step size of the hybrid R algorithm 
and stopping condition for the solver of Kogut- 
Susskind matrix could also affect the value of sus- 
ceptibilities. 

For these reasons the Bielefeld group has con- 
tinued their study pSl, and the JLQCD Collabo- 



ration! 29 1 ^as started their own work last year. 
As one sees in Table |^ run parameters of new 
simulations are chosen to examine the points (i) 
and (ii) above. In order to deal with (iii) Biele- 
feld group worked out the correction formula for 
the case of the single noise vector. They also 
employed the method of multiple noise vectors 
for some of the runs. JLQCD employed the 
method of wall source without gauge fixing [^l| , 
and removed contamination by a correction for- 
mula. At present both groups have accumulated 
(5 — 10) X 10'^ trajectories of unit length with 
a small step size oi 6t = (1 — 1/2) for each 
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Figure 4. Peak height of Xm as a function of 
for spatial sizes L = 8, 12, 16. 

value of [3 and niq. The standard reweighting 
technique is used to find the peak of suscepti- 
bilities. 

In Fig. ^ we plot results of the two groups for 
the peak height of the chiral susceptibility Xm 
for three spatial sizes L = 8, 12 and 16. The 
two results are consistent. A striking feature in 
Fig. ^ is that, except for the largest quark mass 
niq ~ 0.075, the peak height exhibits a signif- 
icant size dependence whose magnitude rapidly 
increases toward small quark masses. In more 
detail, for ruq = 0.0375 and 0.02, the increase is 
largest between L — '& and 12, while the peak 
height is consistent between L = 12 and 16. At 
rUq = 0.01 the peak height increases by a factor 
3 between L = 8 and 16. Data on a 12^ lattice is 
not yet available for this quark mass. 

We should remark that the new data for rUq — 
0.01 for L = 16[^ do not agree with the previ- 
ous results reported in ref. |2^]. For the suscep- 
tibilities Xc and xn for which a direct compari- 
son is possible, new simulations give values which 
are a factor two larger than the old results. It 
is possible that a smaller statistics (2500 trajec- 
tories |2^ as compared to 4400 [|9|), and perhaps 
also a slightly smaller estimate of the critical cou- 
pling (/3c = 5.265 as compared to Pc = 5.266), led 
to an underestimate of susceptibilities in ref. |27) . 

We list the exponents obtained through fits of 
form (6-8) in Table |. While some systematic dif- 
ference appears present for zt^i between the two 
groups, a significant increase of z„i and zt,i from 



Table 3 

Exponents for Nf = 2 
row represents Bielefeld results|28 
those of JLQCDpOl. Results are not yet avail- 



For each exponent first 
and second 



able for entries marked with " 






0{2) 0(4) MF L =8 


L=12 L=16 




0.60 0.54 2/3 0.77(14) 






0.70(10) 


0.63(5) 




0.79 0.79 2/3 0.79(4) 


1.05(8) 0.93(9) 




0.70(3) 


1.01(11) 1.02(7) 


Zt 


0.39 0.33 1/3 






0.65(7) 






0.42(4) 


0.75(12) 0.78(8) 


Zt,a- 


0.63(7) 


0.96(12) 0.86(11) 




0.48(4) 


0.79(14) 0.81(9) 




0.63(7) 


0.94(13) 0.85(12) 




0.47(4) 


0.81(14) 0.82(9) 



L = 8 to L = 12 — 16 is evident, with the val- 
ues for larger sizes sizably deviating from either 
0(4),0(2) or the mean-field predictions. For Zg 
the deviation seems less apparent though full data 
are not yet available. 

A puzzling nature of the values of exponents 
becomes clearer if we translate them into the 
more basic thermal and magnetic exponents yt 
and Uh using the relations. 



- Ml 
Vh 



— , = 1 1 (9) 

Vh Uh Vh 



with c? = 3 the space dimension. The values in 
Table |^ are reasonably consistent with the rela- 



tion Zq -\- Zn 



Zt + 1 which follow from 



We 



observe that Zm ~ 1-0(1) obtained for larger spa- 
tial lattices implies yh ~ 3.0(3) to be compared 
with the 0(4) value 2.49, while yh = d — 3 is ex- 
pected for a first-order phase transition. For the 
thermal exponent we find yt ~ 2.4(3) if we take 
Zt « 0.8(1) or yt « 2.7(3) for Zt « 0.9(1), which is 
substantially larger than the 0(4) value of 1.34. 

One may think of various possibilities for the 
reason leading to these values of exponents. 

(i) The most conventional would be that the 
influence of sub-leading and analytic terms is 
still sizable at the range of quark mass explored. 

(ii) Another possibility, suggested by the value 
yh ~ d for L = 12 and 16, is that a disconti- 
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nuity fixed point with yh = rf ||33| controlling the 
first-order transition along the line = in the 
low-temperature phase is strongly influencing the 
scaling behavior. The transition is of second or- 
der in this case. Whether the deviation of yt from 
any of the expected values can be explained is not 
clear, however, (iii) The transition is of second or- 
der with the exponents close to but not equal to 
d. This would mean a significant departure from 
the universality concepts, stepping even beyond 
the suggestion of mean-field exponents arbitrarily 
close to the critical point |Q. (iv) The transition 
is of first order. In this case, the value of quark 
mass TUq = at which the first-order transition 
terminates would have to be small or even vanish 
since the scaling formula is derived under the as- 
sumption of a transition taking place at a single 
point at = 0. 

Concerning the possibility (iv), results of 
present data examined from finite-size scaling 
point of view are as follows. As we already 
pointed out, Xm for a fixed value of stays 
constant for L = 12 — 16 down to niq — 0.02. 
Results for other susceptibilities exhibit a similar 
behavior. Thus a phase transition does not exist 
for mq_> 0.02 as concluded in the previous stud- 
ies!^ 27 1. At TTiq = 0.01 the susceptibilities in- 
crease by a factor 3 between L = 8 and 16. Runs 
for L ~ 12 are needed to see if the increase is con- 
sistent with a linear behavior in volume expected 
for a first-order transition. 

We have to conclude that scaling analyses of 
susceptibilities carried out so far do not allow a 
definite conclusion. Much further work, possibly 
with a quark mass smaller than has been explored 
so far, is needed to elucidate the nature of the 
chiral transition for Nf = 2. 

3.1.2. Scaling analysis of chiral order pa- 
rameter 

For a second-order transition the singular part 
of the chiral order parameter is expected to scale 
as (qq) = mq^^ 4> (x) where (f>{x) is a function of 
the scaling variable x = {g^^ — g'^^{0))/ml^^^ , 
and l/i5 = 1 - d/yh and l/(3S = yh/vt- A pre- 
vious analysis iQ employing a collection of data 
generated over the years did not find clear sign 
of scaling with the 0(4) exponents. This year 



the MILC Collaboration attempted a more elab- 
orate analysis as part of their study of equation 
of state ||. 

For systems in the 0(4) universality class, the 
scaling function (pix) may be determined in a 
parameterized form through a simulation of the 
0(4) sigma model up to an overall constant for 
the scaling variable a;[^. The result for (j){x) is 
used to fit data for {qq) generated on a 12"^ x 6 
lattice with = 0.025 and 0.0125. Adding an 
analytic term of form mg(co -I- ci/g"^ + €2/3^), 
the fit was found acceptable for 0(4) and also 
for the mean-field scaling function. Extrapolat- 
ing to the limit niq — 0, the results differ signifi- 
cantly between the two cases, however (see Fig. 2 
of ref. H). 

We note that the results of the present MILC 
analysis do not contradict those of susceptibili- 
ties: the quark mass used for this work corre- 
sponds to niq « 0.02 — 0.04 on an S'' x 4 lattice, 
for which case the exponents found from suscepti- 
bilities are similar to the 0(4) values. We further 
remind, however, that the exponents exhibit a sig- 
nificant size dependence. This means that stud- 
ies with larger lattice sizes and smaller rUq are 
required to explore the nature of the two-flavor 
transition from scaling of the chiral order param- 
eter. 

3.2. Restoration of C/a(1) symmetry 

For sufficiently high temperatures topologically 
non-trivial gauge configurations are suppressed, 
leading to restoration of ?/a(1) symmetry. To 
what extent Ua{^) symmetry is restored close to 
the chiral transition is an interesting question. 

Three groups | |35| , |37| , p8[ examined the problem 
using the susceptibility defined by 

XuMD = J {{^{x) ■ m) ~ (5o(x) • 5o(0)))(10) 

which should vanish at mq — if Ua{^) sym- 
metry is restored. In Fig. |5| we plot the rUq de- 
pendence of this quantity obtained by the MILC 
Collaboration!^^ and the Columbia group p7|. 
Both results are taken in the high temperature 
phase corresponding to T/Tc ~ 1.2 — 1.3. While 
the data appear to extrapolate linearly to zero 
at niq = (dotted lines) it is more reason- 
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Figure 5. Susceptibility for C/a(1) symmetry as a 
function of niq/T. 

able to fit with a quadratic dependence X;7a(i) ~ 
a+bmg + 0{m'^) (solid lines) |^ since the suscep- 
tibility is expected to be an analytic function of 
rUq and hence even in ruq in the high temperature 
phase. These fits, which have reasonable lead 
to a non-zero value of XUaW ~ 0- Thus ef- 

fect of anomaly still breaks Ua{^) symmetry just 
above the chiral transition. Similar results were 
reported in ref. p^ . 

The Illinois group calculated the screen- 
ing mass of cr, tt, oq on a 16^ x 8 lattice at 
•niq = 0.00625, employing a noisy estimator for 
the disconnected contribution for the a propaga- 
tor. They found a decrease of the mass splitting 
m,r — across the transition, which, however, 
remains at the level of 20% just above Tc for the 
quark mass employed. They have also shown that 
the disconnected part of the a propagator is dom- 
inated by fermionic modes with small eigenvalues 
induced by instantons. The latter results paral- 
lel those of a previous study carried out for the 
quenched Wilson case||39[]. 

3.3. Energy density and pressure 

The MILC Collaboration completed their 
study of equation of state for the temporal size 
Nt — 6|35| , following their previous work for 
Nt — 4[[40[. They also attempted an extrapola- 
tion of equation of state toward niq — > making 
use of the scaling function computed for the chiral 
order parameter. 

The Bielefeld group fl^ made a measurement of 



energy density for Nf ~ 4 employing the three- 
link improved form of the Kogut-Susskind ac- 
tion]^ together with S(l,2)ti-eG for the gauge ac- 
tion. 

In the results of both groups the energy den- 
sity e/T^ rapidly rises across Tc and stays close 
to the continuum Stefan-Boltzmann value in the 
high temperature phase. In the MILC result a 
bump just above Tc seen for Nt — A almost disap- 
pears for TVt = 6. A similar bump, while observed 
by the Bielefeld group at finite in their Nf = A 
results, is no longer present if terms vanishing in 
the chiral limit is removed. 

These results indicate that e/T"' for full QCD 
is a monotonically increasing function of temper- 
ature similar in shape to that of the pure gauge 
system!^. The behavior of pressure 2>p/T^ is also 
similar, rising smoothly from T k Tc and reach- 
ing e/T^ at T/Tc « 2 - 3. 

4. Phase structure for the Wilson quark 
action 

4.1. Previous phase diagram studies 

A basic concept in the phase structure analysis 
for the Wilson quark action is that of the critical 
line K = Kc{(3) which is usually defined as the 
line of vanishing pion mass. At zero temperature 
this fine runs from w 1/4 at /3 = to = 1/8 
at /3 = oo, and chiral symmetry is expected to 
become restored toward the weak-coupling limit 
along this line[Q. 

At finite temperatures there also exists the line 
of finite-temperature transition K = Kt{P), the 
thermal line. This line starts from the point of 
the deconfinement transition of the pure gauge 
system at K — 0, and moves toward the critical 
line. Naively one would expect the thermal line 
to hit the critical line at some finite f3 = Pet, 
separating the physical region K < Kc{/3) into 
low and high temperature phases. 

Extensive studies have been carried out to ex- 
amine if this expectation is realized p3|-|49||. A 
conceptual issue that arose in the course of stud- 
ies is whether one can naturally define the criti- 
cal line in the high temperature phase since pion 
mass does not vanish in this phase. 

The QCDPAX Collaboration took the vicw||| 
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that the critical hne should be defined by the van- 
ishing of the quark mass iriq at zero temperature, 
where is defined through chiral Ward iden- 
tity |||,|g,|ll. They reported that the crossing 
point Pet with this definition of the critical line 
is located in the region of strong coupling on an 
Nt = 4 lattice, e.g., Pet « 3.9 - 4.0 for Nj = 2. 
For the phase diagram based on this result see 
ref. §. 

This phase diagram, however, has an unsatis- 
factory feature. It has been observed^,^ that 
physical observables do not exhibit any singu- 
lar behavior across the critical line in the high 
temperature phase. This means that the region 
K > Kc{P), usually thought unphysical, is not 
distinct from the high temperature phase, be- 
ing analytically connected to it. Hence one can 
cross from the low- to the high-temperature phase 
through the part of the critical line below /? = Pet , 
which is not a line of finite-temperature transi- 
tion. 

Clearly the phase diagram above does not cap- 
ture the full aspect of the phase structure. Re- 
cent investigations indicate that a more natural 
understanding of the phase structure is provided 
by a different view on the critical line proposed 
by Aoki some time ago|7j. In the following we 
review the phase structure based on this view. 

Let us note that a slightly different phase struc- 
ture has been discussed in ref. |5^. The phase 
structure for general values of Nj up to Nj — 300 
has also been examined recently Q. 

4.2. Spontaneous breakdown of parity- 
flavor symmetry and massless pion 

In order to illustrate the basic idea, let us con- 
sider an effective sigma model for lattice QCD 
with the Wilson quark action with TV/ = 2. The 
effective lagrangian may be written as 

^eff = (V^7f)2 + (V^f7)2 + ^ 5(72 + . . . 

where the coefficients a and b differ reflecting ex- 
plicit breaking of chiral symmetry due to the Wil- 
son term. We know that the pion mass vanishes 
as a = cx Kc—K toward the critical line, while 
a stays massive, i.e., b = > at if « Kc. If 
K increases beyond Kc, the coefficient a becomes 
negative. Hence we expect the pion field to de- 



velop a vacuum expectation value (tt) ^ 0. The 
condensate spontaneously breaks parity and fla- 
vor symmetry. 

Let us note that pion is not the Nambu- 
Goldstone boson of spontaneously broken chiral 
symmetry in this view. Instead it represents 
the massless mode of a parity-flavor breaking 
second-order phase transition which takes place 
aX K = Kc. We expect it to become the Nambu- 
Goldstone boson of chiral symmetry in the contin- 
uum limit, however, as chiral symmetry breaking 
effects disappear in this limit. 

The idea above has been explicitly tested for 
the two-dimensional Gross-Neveu model formu- 
lated with the Wilson action |Q. An analytic so- 
lution in the large N limit shows spontaneous 
breakdown of parity for K > KdP). Another 
important result of the solution is that the criti- 
cal line forms three spikes, which reach the weak- 
coupling limit g = at 1/2K = +2,0,-2. This 
structure arises from the fact that the doublcrs 
at the conventional continuum limit {g, 1/2K) — 
(0, 2) become physical massless modes at 1/2K = 
and -2. 

A close similarity of the Gross-Neveu model 
and QCD regarding the asymptotic freedom and 
chiral symmetry aspects leads one to expect a 
similar phase structure for the case of QCD ex- 
cept that the critical line will form five spikes 
reaching the continuum limit because of differ- 
ence in dimensions 1^]. Evidence supporting such 
a phase structure is summarized in ref. [Q. 

4.3. Finite-temperature phase structure 

For a finite temporal lattice size Nt correspond- 
ing to a finite temperature, the above considera- 
tion can be naturally extend by defining the crit- 
ical line as the line of vanishing pion screening 
mass determined from the pion propagator for 
large spatial separations. 

In Fig. H the critical line for the two- 
dimensional Gross-Neveu model calculated in the 
large N hmit is plotted for Nt = oo,16, 8,4, 2 
starting from the outermost curve and moving to- 
ward inside. The result shows that the location 
of the critical line as defined above depends on 
Nt. Another important point is that the spikes 
formed by the critical line moves away from the 
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Figure 6. Critical line in {g, \/2K) plane for the 
two-dimensional Gross-Neveu model for the tem- 
poral size Nt = oo, 16, 8, 4, 2 (from outside to in- 
side) ||. 

weak-coupling limit as Nt decreases. 

Simulations to examine if lattice QCD has 
a similar structure of the critical line at finite 
temperatures have been made recently for the 
case of iV/ = 211 and 4|5^] on an 8^ x 4 lat- 
tice. The results are summarized as follows: (i) 
For both systems the conventional critical line 
turns back toward strong coupling forming a cusp, 
whose tip is located at /3 ~ 4.0 for Nf — 2 
and /3 « 1.8 for Nf — 4. The cusp repre- 
sents one of five cusps expected for lattice QCD. 
(ii) Parity and flavor symmetry are spontaneously 
broken inside the cusp. Simulations have been 
made for the Nf = 2 system with an exter- 
nal field term SSw — 2KH J2n V'n*75'^3V'ri added 
to the action. Results provide evidence for the 
behavior lim.H^o{''Pl5'''34') 7^ of the parity- 
fiavor order parameter and vanishing of tt^ mass 
MmH^o ™7r± = expected inside the cusp|]5^. 

Concerning the relation between the thermal 
line and the critical line, we recall that the pion 
mass vanishes all along the critical line. This 
suggests that the region close to the critical line 
is in the cold phase even after the critical line 
turns back toward strong coupling, and hence the 
thermal line cannot cross the critical line. Since 
numerical estimates show that the thermal line 
comes close to the turning point of the cusp, the 
natural possibility is that the thermal line runs 
past the tip of the cusp and continues toward 
larger values of K. Results of measurement of 
thermodynamic quantities provide support of this 



Figure 7. Phase diagram for Nf ~ 2,3,4 on an 
TVf = 4 lattice. 

view[^, although the possibility that the ther- 
mal line touches the critical line at a point ||4^ 
cannot be excluded. 

In Fig. we summarize presently available re- 
sults for the phase structure on an iVi = 4 lattice 
for Nf = 2, 3 and 4. Solid lines represent the 
critical line estimated from the pion mass. For 
the case of Nf — 3 results confirming the cusp 
structure arc not yet available. Open and solid 
squares show simulation results for the location 
of the thermal line. Dotted lines are smooth in- 
terpolation, extended beyond the cusp following 
the discussion of the previous paragraph. The 
differe nce of open and solid squares is discussed 
in Sec. 4.4 below. 

Let us remark that the conventional zero tem- 
perature critical line runs close to the lower part 
of the finite-temperature critical line in Fig. |^ and 
continues toward weak coupling. The thermal 
line therefore has to cross the zero-temperature 
critical line. This represents the crossing point 
f3ct reported by the QCDPAX Collaboration[||. 

We emphasize, however, that the zero- 
temperature critical line does not represent a line 
of singularity of the finite-temperature partition 
function. This naturally explains the absence of 
singular behavior of observables across the zero- 
temperature critical line mentioned in Sec. 4.1. 



4.4. Nf dependence of order of transition 

The sigma model analysis in the continuum|^ 
suggests that the chiral phase transition is of sec- 
ond order for Nf — 2 and of first order for Nf > 3. 
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Indeed strong first-order signals have been ob- 
served for the case oi Nf = and 4|56| away 
from the critical line, as shown by solid squares 
in Fig. 1^, in contrast to a crossover behavior rep- 
resented by open squares seen for Nf — 2. How- 
ever, the first-order transition for Nf = 4 weak- 
ens closer to the critical line, apparently turning 
into a smooth crossover before reaching the region 
around the cusp of the critical line as indicated 
by open rectangles |Q. While parallel data are 
not yet available for Nf = 3, results of the QCD- 
PAX Collaboration!^^ also appears to indicate a 
weakening of the first-order transition. 

A possible reason for this unexpected behav- 
ior is that breaking of chiral symmetry due to 
the Wilson term, which becomes stronger as (3 
decreases along the thermal line, smoothens the 
first-order transition. Another possibility is that 
the first-order transition for Nf — 3 and 4 ob- 
served so far is a lattice artifact sharing its origin 
with the sharpening of the crossover at /3 « 5.0 
found by the MILC Collaboration for A^/ = 2@. 
Some support for this interpretation is given by 
a recent study of the QCDPAX Collaboration for 
the Nf = 3 system with an improved gauge ac- 
tion!^. So far they have not found clear first- 
order signals in the region where the plaquette 
action shows a clear first-order behavior. 

In either case, if chiral transition in the con- 
tinuum is indeed of first order for A^^ — 3 and 
4, it will emerge only when the cusp moves suffi- 
ciently toward weak-coupling with an increase of 
the temporal size Nt- 

4.5. Continuum limit 

We expect the cusp of the finite-temperature 
critical line to grow toward weak coupling as Nt 
increases. In the limit A^t = oo it should con- 
verge to the zero-temperature critical line which 
reaches /3 = oo. Since the thermal line is located 
on the weak-coupling side of the cusp for a finite 
Nt, it will be pinched by the tip of the cusp at 
K) = (oo, 1/8) as Nt oo. We expect chiral 
phase transition in the continuum to emerge in 
this limit. In order to extract continuum proper- 
ties of the chiral transition, we then need a sys- 
tematic study of thermodynamic quantities in the 
neighborhood of the thermal line when it runs 



close to the tip of the cusp as a function of Nt- 

Simulations, however, indicate that the cusp 
moves only very slowly as Nt increases. For the 
Nf = 2 case, current estimates of the position 
of the tip of the cusp is /? w 4.0 for Nt = 4 [^ 
4.0 - 4.2 for Nt = 6|§, 4.2 - 4.3 for A^f = 8|56| 
and 4.5 — 5.0 even for A't = 18|^^. A recent work 
also reports an absence of parity-broken phase 
above (3 — 5.0 on symmetric lattices up to the 
size lO^'Q. For Nf = A the values are even 
lower:/? « 1.80 for A^f = 4 and 2.2 - 2.3 for 
Nt — 8[|6|. These estimates indicate that a very 
large temporal size will be needed for the cusp 
to move into the scaling region {e.g., (3 > 5.5 for 
Nf = 2) as long as one employs the Wilson quark 
action together with the plaquette action for the 
gauge part. 

We emphasize that this result has an impor- 
tant implication also for spectrum calculations at 
zero temperature. Since the location of the cusp 
is determined by the smaller of the spatial and 
temporal size, the critical line will be shifted or 
may even be absent unless lattice size is taken 
sufficiently large. Therefore hadron masses cal- 
culated on a lattice of small spatial size and ex- 
trapolated toward the position of the critical line 
might involve significant systematic uncertainties. 

4.6. Studies with improved actions 

The problems discussed above indicate the 
presence of sizable cutoff effects when the Wil- 
son quark action is used in conjunction with the 
plaquette action. A way to alleviate this problem 
is to employ improved actions. This approach 
has been pursued by the QCDPAX Collabora- 
tion |8||5^, replacing the plaquette action with an 
improved gauge action RG(1,2) |l^ . This year the 
MILC Collaboration reported simulations with 
the action SUWtadpoie for the gauge part and 
the tadpole-improved clover action for the quark 
part[Q. Results with the tree- level clover ac- 
tion keeping the plaquette action are also avail- 
able!^^. Thus there are data for four types of 
action combinations, unimproved and improved 
both for the gauge and quark actions, to make a 
comparative study of improvement. 

An indication from such a comparison is that 
improving the gauge action substantially removes 
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cutoff effects. An inflection of tlie critical line seen 
for the plaquette action at /3 w 4 — 5 becomes 
absent with improvement of the gauge action js), 
while it still seems to remains if only the Wilson 
quark action is replaced by the clover action 
Also an intermediate sharpening of the thermal 
transition seen for the plaquette action at /3 « 
5.0 is not observed for improved actions . 

Another point to note is that the lattice spacing 
at the coupling constant where the thermal line 
approaches the critical line has a similar value 
TOpfl Ri 1 on an. Nt = 4 lattice for all of the four 
action combinations. This means that studies of 
physical quantities are needed to assess reduction 
of cutoff effects with improved actions. Interest- 
ing results have already been obtained for scaling 
of the chiral order parameter |^,^, and work with 
the critical temperature is being pursued pjjs^]. 

5. Results in finite density studies 

It has long been known that the quenched ap- 
proximation breaks down for a non-zero quark 
chemical potential /x in that a transition takes 
place at /X « m7r/2 rather than at /i « mAr/3]6^, 
|6^ . While the importance of the phase of 
the quark determinant has been made clear, 
the mechanism how the quenched approximation 
breaks has not been fully explained. 

Recently Stephanovp^, employing a random 
matrix model of the quark determinant []65| and 
a replica formulation of quenched approximation, 
traced back the failure of the quenched approxi- 
mation to the non-uniformity of the limit of the 
replica number n —^ for /i = and /i > 0. He 
has also shown that the quenched approximation 
is valid for the theory in which a quark x in the 
conjugate representation is added to each quark 
q. Formation of a condensate (xq) having a unit 
baryon number and a mass m^q « in such 
a theory explains the occurance of transition at 

Barbour and collaborators reported new results 
in full QCD simulations^^. With the method 
of fugacity expansion runs were carried out 
for four flavors of quarks on 6^ and S'' lattices at 
/3 = 5.1 with the Kogut-Susskind quark action. 
They found an onset of non-zero baryon number 



at a small value of /x, e.j., \ic ~ 0.1 at mq — 0.01. 
For comparison the MTc Collaboration reported 
TOjv = 1.10(6) and = 0.290(6) at a slightly 
larger couphng of /3 = 5.15 at rriq = 0.01[|68|. 

It is not yet clear if these results mean that an 
early onset of transition fic ~ also holds 

for full QCD or reflect computational problems 
of the method employed for the simulation. 

6. Conclusions 

Much work has been made in finite tempera- 
ture studies of lattice QCD encompassing a num- 
ber of subjects during the last year. 

Tests of improved actions made for the pure 
gauge system indicate a possibility that accurate 
results for thermodynamics in the continuum may 
be obtained with simulations carried out with a 
moderately large temporal size. 

In full QCD studies much progress has been 
made in understanding the phase structure for 
the Wilson quark action. On the other hand, 
new problems have also been encountered, mak- 
ing it necessary to reexamine conclusions reached 
in previous studies. These are the unexpected val- 
ues of exponents for Nf — 2 found in scaling stud- 
ies of susceptibilities with the Kogut-Susskind 
quark action, and the flavor dependence of or- 
der of chiral transition with the Wilson quark ac- 
tion. Elucidating these problems is important for 
reaching an understanding of the nature of chiral 
phase transition, which is consistent between the 
Kogut-Susskind and Wilson quark actions. 

Some progress has been made in QCD at finite 
density. A puzzling result reported from the lat- 
est simulation shows, however, that we are still 
far from understanding this difficult subject. 

In closing we point out that most work in full 
QCD during the past several years have concen- 
trated on the case of Nf degenerate quarks, espe- 
cially for Nf = 2. While a variety of basic prob- 
lems we have encountered for this case has to be 
clarified with further work, we should also recall 
that nature corresponds to the case of Nf ~ 2 + 1 
with a heavier strange quark. A delicate change 
of phase that might possibly result from its pres- 
ence, as suggested in the continuum sigma model 
analysis 1^, makes it important to enlarge previ- 
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ous studies [ p7||6S| Ji^ into a systematic effort in 
this direction. 

Acknowledgements 

I would like to thank S. Aoki, F. Beinlich, 
T. Blum, W. Bock, G. Boyd, M. Creutz, N. 
Christ, T. DcGrand, C. DeTar, M. Fukugita, Y. 
Iwasaki, K. Kanaya, T. Kaneko, F. Karsch, E. 
Laermann, P. Mackenzie, M. Okawa, D. Tous- 
saint, M. Wingate and Y. Yoshie for communicat- 
ing their results and for discussions. I would also 
like to thank S. Aoki, M. Fukugita, Y. Iwasaki, 
K. Kanaya and M. Okawa for comments on the 
manuscript. This work is supported in part by 
the Grant-in- Aid of the Ministry of Education, 
Science and Cuhure (Nos. 04NP0801, 08640349). 

REFERENCES 

1. For a review and references, see, A. Ukawa, 
Lattice 89, Nucl. Phys. B (Proc. Suppl.) 17 

(1990) 118. 

2. For a representative work, see, G. Boyd, J. 
Engels, F. Karsch, E. Laermann, C. Legeland, 
M. Liitgemeier and B. Petersson, Nucl. Phys. 
B469 (1996) 419. 

3. For a summary as of 1990, see, S. Gottlieb, 
Lattice 90, Nucl. Phys. B (Proc. Suppl.) 20 

(1991) 247. 

4. F. Karsch, Phys. Rev. D49 (1993) 3791; F. 
Karsch and E. Laermann, Phys. Rev. D50 
(1994) 6954. 

5. R. D. Pisarski and F. Wilczek, Phys. Rev. 
D29 (1984) 338; F. Wilczek, Int. J. Mod. 
Phys. A7 (1992) 3911; K. Rajagopal and F. 
Wilczek, Nucl. Phys. B399 (1993) 395. 

6. For previous reviews, see, Y. Iwasaki, Lattice 
94, Nucl. Phys. B (Proc. Suppl.) 42 (1995) 96; 
K. Kanaya, Lattice 95, Nucl. Phys. B (Proc. 
Suppl.) 47 (1996) 144. 

7. S. Aoki, Phys. Rev. D30 (1984) 2653; Phys. 
Rev. Lett. 57 (1986) 3136; Nucl. Phys. B314 
(1989) 79. 

8. Y. Iwasaki, K. Kanaya, K. Kaya, S. Sakai and 
T. Yoshie, Lattice 94, Nucl. Phys. B (Proc. 
Suppl.) 42 (1995) 502; Lattice 95, Nucl. Phys. 
B (Proc. Suppl.) 47 (1996) 515. 



9. T. DeGrand, A. Hasenfratz, P. Hasenfratz 
and F. Niedermeyer, Nucl. Phys. B454 (1995) 
587, 615. 

10. M. Blatter and F Niedermeyer, BUTP-96/12 
( [hep-lat/96050T7| ). 

11. A. Papa, BUTP-96/13 ( |hep-lat/9605004| ). 

12. Y. Iwasaki et aZ.(presented by T. Kaneko), 
this volume, UTHEP-3 41 (|hep-lat/9608090D 

13. W. Bock, this volume ( |hep-lat/9608103D ; W. 
Bock, Y. Iwasaki, K. Kanaya and Yoshie, in 
preparation. 

14. G. CeUa, G. Curci, A. Vicere and B. Vigna, 
Phys. Lett. B333 (1994) 457. 

15. F. Karsch, t his volume, BI-TP 96/34 (^ 



lat/9608047) 



16. F. Beinlich, F. Karsch and E. Laermann, 
Nucl. Phys. B462 (1996) 415. 

17. D. W. Bhss, K. Hornbostel and G. P. Lepage, 



18 



SMUHEP 96-05 ( |hep-lat/960504lD . 
Y. Iwasaki, Nucl. Phys. B258 (1985) 141; 
Tsukuba preprint UTHEP-118 (1983), un- 
published. 

19. G. P. Lepage and P. Mackenzie, Phys. Rev. 
D48 (1993) 2250. 

20. M. Liischer and P. Weisz, Phys. Lett. B158 
(1985) 250 and references therein. 

21. M. Alford, W. Dimm, G. P. Lepage, G. 
Hockney and P. Mackenzie, Phys. Lett. B361 
(1995) 87. 

22. S. Collins, R. Ed wards, U. Heller and J. 
Sloan, this volume ( |hep-lat/960802l| ). 

23. T. DeGrand, private communication. 

24. P. de Forcrand, G. Schierholz, H. Schneider 
and M. Teper, Phys. Lett. B160 (1985) 137. 

25. Y. Iwasaki et al, Phys. Rev. D46 (1992) 
4657;Phys. Rev. D49 (1994) 3540;Y. Aoki and 
K.Kanaya, Phys. Rev. D50 (1994) 6921. 

26. M. Fukugita, H. Mino, M. Okawa and A. 
Ukawa, Phys. Rev. Lett. 65 (1990) 816; Phys. 
Rev. D42 (1990) 2936. 

27. F. R. Brown, F. P. Butler, H. Chen, N. H. 
Christ, Z. -H. Dong, W. Schaffer, L. I. Unger 
and A. Vaccarino, Phys. Rev. Lett. 65 (1990) 
2491; A. Vaccarino, Lattice 90, Nucl. Phys. B 
(Proc. Suppl.) 20 (1991) 263. 

28. G. Boyd, F. Karsch, E. Laermann and M. 



Oevers, BLTP 96/27 (|hep-lat/9607046D ; in 
progress (private communication from F. 



14 



Karsch and E. Laermann). 

29. JLQCD Collaboration, S. Aoki et ai, in 
progress. 

30. A. Kocic and J. B. Kogut, Phys. Rev. Lett. 
74 (1995) 3109; Nucl. Phys. B455 (1995) 229. 

31. Y. Kuramashi et ai, Phys. Rev. Lett. 72 
(1994) 3448. 

32. I. R. Mcdonald and K. Smger, Disc. Faraday 
Soc. 43 (1967) 40; A. M. Ferrenberg and R. 
Swendsen, Phys. Rev. Lett. 61 (1988) 2635. 

33. B. Nienhuis and M. Nauenberg, Phys. Rev. 
Lett. 35 (1975) 477. 

34. C. DeTar, Lattice 94, Nucl. Phys. B (Proc. 
Suppl.) 42 (1995) 73. 

35. MILC Collaboration (presented by C. Detar 
and T. Blum), this volume ( |hep-lat/9608026D . 

36. D. Toussaint, AZPH-TH/96-13 (hep-lat 
/ 9607084) 

37. N. Christ, this volume. 

38. J. Kogut, J.-F. Lagac and D. Sinclair, 
this volume, ANL-HEP-CP-96-52 (hep- 
lat/9608128). 

39. S. Itoh, Y. Iwasaki and Y. Yoshie, Phys. Rev. 
D36 (1987) 527. 

40. T. Blum, S. Gottlieb, L. Karkkiiinen and D. 
Toussaint, Phys. Rev. D51 (1995) 5153. 

41. S. Naik, Nucl. Phys. B316 (1989) 238. 

42. M. Bochicchio, L. Maiani, G. Martinelh, G. 
Rossi and M. Testa, Nucl. Phys. B262 (1985) 
331. 

43. M. Fukugita, S. Ohta and A. Ukawa, Phys. 
Rev. Lett. 57 (1986) 1974. 

44. A. Ukawa, Lattice 88, Nucl. Phys. B (Proc. 
Suppl.) 9 (1989) 463. 

45. R. Gupta, A. Patel, C. Baillie, G. Guralnik, 
G. Kilcup and S. Sharpe, Phys. Rev. D40 
(1989) 2072. 

46. K. M. Bitar, A. D. Kennedy and P. Rossi, 
Phys. Lett. B234 (1990) 333. 

47. HEMCGC CoUaboration, K. M. Bitar et at, 
Phys. Rev. D43 (1991) 2396. 

48. Y. Iwasaki et ai, Phys. Rev. Lett. 67 
(1991) 1491; Ph ys. Rev. Lett. 69 (1992) 
21; UTHEP-300(ihep-lat/9504019|) ;UTHEP- 
304 (|hep-lat/9505017|) . A Ml account is re- 
ported in, UTHEP-334 (|hep-lat /9605030|) . 

49. C. Bernard, M. Ogilvie, T. A. DeGrand, C. 
DeTar, S. Gottlieb, A. Krasnitz, R. L. Sugar 



and D. Toussaint, Phys. Rev. D46 (1992) 
4741; C. Bernard et aLPhys. Rev. D49 (1994) 
3574; Lattice 93, Nucl. Phys. B (Proc. Suppl.) 
34 (1994) 324; T. Blum et ai, Phys. Rev. 
D50 (1994) 3377; C. Bernard et ai, Lattice 
94, Nucl. Phys. B (Proc. Suppl.) 42 (1995) 
451. 

50. S. Itoh, Y. Iwasaki, Y. Oyanagi and T. Yoshie, 
Nucl. Phys. B274 (1986) 33. 

51. L. Maiani and G. Martinelh, Phys. Lett. B178 
(1986) 265. 



55 



56 



52. M. Creutz, |hep-lat/9608024 . 

53. Y. Iwasaki et ai, (presented by K. Kanaya 
and S. Kaya), this volume (hep/lat-9608125). 

54. S. Aoki, Prog. Theor. Phys. 122 (1996) 179 
and references cited therein. 
S. Aoki, A. Ukawa and T. Umemura,Phys. 
Rev. Lett. 76 (1996) 873; Lattice 95, Nucl. 
Phys. B (Proc. Suppl.) 47 (1996) 511. 
S. Aoki, T, Kaneda, A. Ukawa and T. 
Umemura, this volume. 

57. Y. Iwasaki et ai, in progress. 

58. MILC Collaboration ( presented by M. 
Wingate), this volume ( |hep-lat/9607085| ). 

59. S. Aoki, U. Ukawa and T. Umemura, unpub- 
lished (1995).^ 

60. K. M. Bitar, |hep-lat/9602027 . 

61. Y. Iwasaki et ai, UTHEP-344 (hep- 
lat /9609022). 

62. I. Barbour et ai, Nucl. Phys. B275[FS17] 
(1986) 296. 

63. M. -P. Lombardo, J. B. Kogut and D. Sinclair, 
Phys. Rev. D54 (1996) 2303 and references 
cited therein. 

64. M. A. Stephanov, Phys. Rev. Lett. 76 (1996) 
4472; this volume ( |hep-lat/960706C| ). 

65. E. V. Shuryak and J. J. M. Verbaarschot, 
Nucl. Phys. A560 (1993) 306. 

66. I. M. Barbour, J. B. Kogut and S. E. Morrison 
(presented by S. E. Morrison), this volume 
( |hep-lat/9608057D . 

67. I. M. Barbour and Z. Sabeur, Nucl. Phys. 
B342 (1990) 269. 

68. MTc Collaboration, R. Altmeyer et ai, Nucl. 
Phys. B389 (1993) 445. 

69. J. B. Kogut and D. K. Sinclair, Phys. Lett. 
B263 (1991) 101. 



